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Abstract We present a general framework for the model-independent decomposition of the fully differential 
cross section of the reactions e + e _ — > 7* — > 7r°7r°7 and e + e~ — > 7* — > ir a rff, which can provide important 
information on the properties of scalar mesons: /o(600), /o(980) and ao(980). For the model-dependent 
ingredients in the differential cross section, an approach is developed, which relies on Resonance Chiral The- 
ory with vector and scalar mesons. Numerical results are compared to data. The framework is convenient 
for development of a Monte Carlo generator and can also be applied to the reaction e + e~ — > 7* — s> 7r + 7r~7. 



1 Introduction 



Despite extensive studies during last decades, physics of 
the light scalar mesons a (980) (I G {J PC ) = l-(0++)), 
/o(980) and / (600) = a (I G (J PC ) = 0+(0++)) is far from 
complete understanding. In particular, there are doubts 
whether simple quark model can explain their properties, 
see, e.g., the review in pQ. 

The dominant decay channels of scalar mesons are 
known to be ir + ir~, 7r°7r° for the /o(980) and a me- 
son, and 7r°?7 for the ao(980) meson. Much experimen- 
tal attention has already been paid to the radiative de- 
cays of the (f> meson: 0(1020) — > r yaQ — > -77777 and 
0(1020) -> 7/0 (or 7 cr) -> 77T7T [US] (see also the KLOE 
summary in [5] and results from Novosibirsk (TUHMH])- Such 
measurements are a good source of information about the 
scalar meson properties [10]. Various models have been 
proposed to describe these decays, [TrJl[TTl[T^lll3l[Tl] . to 
mention a few. The calculated decay widths turn out to 
be very sensitive to model ingredients, however, the ex- 
perimental data is still insufficient to unambiguously dis- 
criminate between the models. 

In the case of the neutral final state (FS), i.e., 7r 7r°7 
and 7r°777, the cross section is determined solely by final- 
state radiation (FSR) mechanism, since there is no initial- 
state radiation (ISR) contribution resulting in the same 
final state. Despite the lower value of the cross section, 
compared to the charged pion case (e + e _ — > tt + tt^j), pro- 
cesses with the neutral-meson FS are an invaluable source 
of information on complicated hadron dynamics. 



In this paper we describe the differential cross section 
of the e + e~ annihilation to a pair of neutral pscudoscalar 
mesons and one photon in the FS, 

e+( P+ ) e -(p_) 7* PiiPi) ft(pa) 7(A). (1) 

The pseudoscalar mesons (J PC = '") are denoted by 
P\P-2 = 7r°7r° and ir°rj. In Section [5] we present a formalism 
for a differential cross section, which is the main task of 
this paper. We provide more general formulae in compar- 
ison with Refs. [15|ll6lH7] . namely, the non-integrated ex- 
pressions are given as well as those integrated over the an- 
gles. It gives a convenient ground to implement the results 
in the Monte Carlo generators, e.g., in FASTERD [18] 
(based on the general structure given in Ref. [IS]) or 
PHOKHARA [TO]. 

Our framework is consistent with symmetries of the 
strong and electromagnetic interactions. It incorporates a 
model-dependent description of the FSR only through the 
explicit form of the Lorentz-invariant functions /i,2,3 and 
has a model-independent tensor decomposition. 

In Sections [3] and 0] we calculate the FS hadronic ten- 
sor. It is the second goal of the paper to provide such a 
description in terms of functions /i,2,3- Our model relies on 
the Lagrangian of Resonance Chiral Theory (R%T ) J^U] ■ 
The R%T is a consistent extension of Chiral Perturba- 
tion Theory to the region of energies near 1 GeV, which 
introduces the explicit resonance fields and exploits the 
idea of resonance saturation. One of the advantages of the 
R^T Lagrangian at leading order (LO), which makes it 
convenient for the present study, is that, having a good 
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predictive power, it contains very few free parameters 
compared with other phenomenological models. In order 
to get good agreement with data, we release a rigor of 
RxT and include some SU(3) symmetry breaking effects 
(e.g., use realistic masses of vector mesons) and mixing 
phenomena (e.g., a G-parity-violating 4>unr° transition). 

The loop contributions follow from the model La- 
grangian. For example, the kaon loop in the 0/o7 
transtion, which is often considered as a pure phenomenol- 
ogy manifestation, in the present model is a direct conse- 
quence of the RxT Lagrangian. In order to simplify the 
formulae, some numerically irrelevant loop contributions 
are omitted. In addition, the resonance exchanges in the 
loops are not considered to avoid problems with renormal- 
izability. 

We consider in detail the following intermediate states 
with scalar and vector resonances, which lead to the same 
FS FiP 27 : 



scalar decay, (Section [3]) 
7* -> S 7 -> P1P27 (2) 
7* -> V -> 57^ PiP 2 7 

vector contribution, (Section 2]) 
7* -> VP h2 -> P1P2I (3) 

7 * -> V a -> V6Pl, 2 -> AP27 



e + e 



e + e 



e + e 




Figure 1. Generic scheme for electron-positron annihilation 
into two particles with final state radiation 



The cross section of the FSR process can be written as 
1 



dap = 



2s(2tt) ; 



■C 



12 



x / 5 4 {Q- Pl - P2 -k)\M FSR \ 



; d 3 p 1 d 3 p 2 d 3 k 



= C 12 N 
1 



J \M F srs\ 2 d cos 6 depdm^ dp 



&E1E2 

(5) 



TV = 



(2tt) 4 64s 2 ' 

where s = Q 2 , 9 is the azimuthal angle, <f> is the polar angle 
of the photon and m\ 1 = (k+pi) 2 . The factor C12 = 1/2 
for 7r°7r° in the final state and C\ 2 = 1 for ir°r]. The matrix 
element Mfsr is 



71/ 



FSR 



(6) 



where S ( J PC = ++ ) is an intermediate scalar meson 
(S = fo, o for ttqttq FS and S = clq for 7Tor;). Only the 
lowest nonet of vector mesons (V, V a , Vb — p, uj and </)) 
is taken into account. 

We are interested in the center-of-mass energy ^/s 
range from the threshold up to M$. This framework may 
also be used in a somewhat dedicated case of y 7 ! = M^, 
giving, e.g., the 4> radiative decay description. 

For the quantitative illustration of our approach, in 
Section [5] we show the numerical results for the values 
of -^/s = 1 GeV and -^/s = M$. The meson-pair invari- 
ant mass distributions are of interest, and for y/s = 
they are compared with available results from KLOE. We 
demonstrate the interplay of the contributions © and 
Conclusions follow in Section [6] 



2 General structure of the FSR cross section 



For a generic reaction e + e — s- 7P1P2 we define 4- 
momenta as shown in Fig. [T] 



p = pi+p 2 , l=Pl~P2, 
Q = P+ + P- = k + pi + p 2 . 



(4) 



The masses of pseudoscalars are m(Pi) = mi, m(P 2 ) = 
m 2 . 



where e = ^J^kol « y/Air/137 w 0.303 and the FSR ten- 
sor M^ v can be decomposed into three gauge-invariant 
independent tensors: 



(7) 



uv 



uv 



WQ V - g^k ■ Q, 

k ■ l{l»Q v - g^Q ■ I) + l v (WQ -l-Fk- Q), 
Q 2 {g^k ■ I - m v ) + Q^(l v k -Q~Q u k-l) 



with the Lorcntz-invariant functions 

/, = fi(Q 2 ,k-Q,k-l), 



(8) 



i = 1, 2, 3. If mi = ni2, these tensors coincide with those of 
Rcf . [TM2T[ . One may also find a similar approach in jTTJ 
22,23]- We emphasize that the decomposition (J7|) is model 
independent; the model dependence is contained in an ex- 
plicit form of functions /j only. Notice that the scalar 
products can be written in terms of the invariant masses: 



k-Q = (s-p 2 )/2, 
k ■ I = vn\ 1 — m\ - 
Q-l = k-l + sS/2, 



k-Q, 



(9) 



where S = 2{m\ — m 2 )/s. 



For the matrix element squared and averaged over the 



e + e polarizations we obtain 

„6 



|M 



FSR\ 



a 11 |/ 1 | 2 + 2a 12 Re(/ 1 / 2 *)+a22|/2| 



2 a 23 Re(hfa) + a 3 3 |/ 3 | 2 + 2a 13 Pe(/ 1 /*) 



(10) 
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aik = (xflw - P+pP- P - p+pP-^r^T^g^x, (11) 



equal to 



of the final-meson phase space we have 

da _ a 3 C r2 
dm^dp 2 32s 



(A 11 |/i| 2 + 2A 12 Re(/ 1 / 2 *) + A 22 |/ 2 | 



+2A 23 i?e(/ 2 / 3 *) + ,4 33 |/ 3 | 2 + 2A 13 i?e(/ 1 / 3 *)) , (14) 



a 22 = I s^ih + i 2 ) 2 + 4l 2 ( Ul 2 (s 2 + s(ti + t 2 ) + 4) 
8 

+U 2 2 (s 2 + S(tl + t 2 ) + t\) 

+2 Ul u 2 (s 2 + s(h + t 2 ) - ht 2 )) 
+8s(u 2 + u 2 2 ){ui +u 2 ) 2 

- (4u 2 1 +4:ul + l 2 (2s + t 1 +t2)) : 
s 3 S 2 

+ (l 2 s + 2ul + 2ul)—-, 



s 2 {ui + u 2 )S 



a.33 = — ^{hhl 2 + 2(mi + u 2 )(u 2 ti + u\h) 
-5s(u 2 ti + uit 2 )), 



where 








A n = 


Ax 
~3~' 






A 12 = 




X 2 ) 2 +X 2 (<7 - 


1 + x) - 25(xi - x 2 ) + S 2 


A u = 




-X2-5) 




A 23 = 


2s 2 
— (si 


— X2)(S — X{ - 




A22 = 


s 2 

Vk»- 

+2x 2 (a- 


x 2 ) 4 + 2( Xl - 
- 1 + xf 


x 2 ) 2 {l - x)(a - l + x) 




-2S(xi - 


X2) ({xi - x 2 


f + (a - l + x)(xi +x 2 )) 




+S 2 ((xi 


-x 2 ) 2 + 2{a- 


-l + x))}, 


A33 = 


2s 2 

-3-[(*i- 


-x 2 ) 2 {l + x)- 


- x 2 {a - l + x) 



1 

«12 = 77 



sl 2 (h + i 2 ) 2 + 4u 2 (s 2 + st 2 + t 2 2 ) 
+4u 2 ,(s 2 + sh + < 2 ) + 4uiu 2 (2s 2 + s(ti + t 2 ) - 2tit 2 ) 
+2s 2 (tiu 2 + t 2 ui + 2s(m + u 2 )) S + s 4 S 2 



+S(S-(2 + x){x 1 -x 2 ))], 



(15) 



and 



Ol3 



023 



(ui + u 2 )(sti + st 2 + tits) - u\t\ - u 2 t 2 



X2 



s — p 
s 

2E 2 

^ = 



2Ex _ p 2 + m 2 7 - m 2 . 



s + m 2 — m 2 



G = 2(m 2 + m 2 ) 



For the case mi = m 2 Eq. (fl~2"j) reduces to Eq. (17) 
of Ref. [15]. Also the results ([H)). (fT5j) coincide with 
Eqs. (2.7), (2.8) of [in]. However, for an MC generator, 
the expressions ([5]) and (|T0]) with coefficients aj& are more 
convenient than (HU). 

+<5s («iu 2 (4s + ti + i 2 ) + u\(2s - t 2 ) + u 2 (2s - t 2 )) Integrating Eq. ([II)) over m 2 one obtains the distri- 



/ 2 (uit 2 - w 2 ti)(ti - t 2 ) - 2s( Ul + u 2 f 



+2{m + u 2 ){ui - u 2 )(t 2 ui ~ u 2 t\) 



8 2 



where 



3 K 


+ u 2 ) 




h = 


(P-- 


k) 2 - m 2 = 


h = 


(p+- 


k) 2 -m 2 e = 


Ul = 


l-p- 


u 2 = l-p+ 



(12) 



bution of the invariant mass \Jp 2 of two pscudoscalar 



-2j>- • fc, 



da 
d\i p 2 



2 VP 



dm 



1- 



(ZfT 



dmij dp 2 



(17) 



o_i_ • fc, The bounds of integration over m 2 7 at the fixed value of 

(13) P 2 are determined by 

For numerical calculations the relation I 2 — 2(m 2 +m|) — r 2 \ _ s {p 2(T + s $) 

p z may be useful. ' ' 

The Eqs. ([5J and (|TD|) . with the explicit expres- 
sions ([12]) and (fT5|) , fix the whole model- independent part 

of the differential cross section. It is worth illustrating a 2 
relation of these formulae to the partial differential cross 

section. Taking into account the corresponding factors and At the </>-meson peak (s = M?) one can present the 
integrating the coefficients aik over the angular variables results in terms of the branching ratio for the cf> — > P\P 2 ^ 



Ap 2 



s — p 



/ sa s 2 S 2 \ , 

1± V 1 "7 + v i (18) 
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Figure 2. Scheme of e + e — > S7 — > P1P2J subprocess 

decay, which is related to the cross section as follows: 
dB{<j> -> P1F27) _ M l 



dy P 2 



12ttB(0 -> e+e") 
da{e+e- -> P1P27) 

^ 7^ 5 



(19) 



where the 4> — > e + e~ branching ratio B{<f> — > e + e _ ) is 
used. In the context of this paper, a calculation of this 
branching ratio is useful for comparison of model predic- 
tions with available data. 



3 Scalar contribution 

In this Section we consider in detail the transition ampli- 
tudes 



7 
7' 

7* -> a 7 



/07 -> 7r°7r°7, 



(77 — > 7r°7r°7, 

7T°777 



(20) 

for the 7r°7r°7 and 7r°?77 final states, respectively. They 
contibute to e + e~ — > S7 — > -P1P27 as illustrated in Fig. [21 
To describe the processes ((2"0"j) we use the Lagrangian of 
R-XT [20] at the linear-in-resonance level, following [TJl 
[24] . The basic features of the Lagrangian framework of 
the RxT are sketched in | Appendix A| We emphasize that 
both light isoscalar scalar resonances, /o and a are in- 
cluded in the formalism in a natural way. Throughout this 
section we work in the tensor representation for spin-1 par- 
ticles [2"0"U2"5] . In the present work we take into account the 
pscudoscalar decay constants splitting (/ w 7^ fx) which 
was discussed in the same context in Ref. [24] , 

The interaction of pseudoscalars with the photon field 
in RxT is identical to the scalar QED. We shall now 
discuss the interaction terms of the Lagrangian (|A3[) rele- 
vant to the processes ([20j) (cf. [H]). For the vector mesons 
in the even-intrinsic-parity sector one has 



£~ f v = eF v F^(±p% + ^ 



3V2 



(21) 



Cvpp = iG v [— (2 p^d^+d^-) 

JK 

+ - V2^)(d»K°d»K )] , (22) 

JK 



J 7T 



■^■^B" {pl v + cu^ - V2^ v ) K+K- 
2eG v 



B v p\ v (tt+^tt- + 7r-^7r+) 



J IT 

eGy jji, 

f 2 
■Ik 

x {K+d^K- + K-d^ l K+) , 



B" (plv+u^-y/2^) 



(23) 



where F al3 stands for the electromagnetic field tensor and 
V^ v for the vector field in the tensor representation, Fy 
and Gv are the model parameters (see |Appendix B| for 
numerical values). Vertex functions for Eas. (|21[) ~(|23 p are 
shown in Table [TJ 

The Lagrangian terms for scalar and pscudoscalar me- 
son interactions, which follow from (|A5[1 are 



Cs 



s 



1 gSirir -> 2 . 1 3Sr,r, 2 , 1 
Jtx Z Jit z Jtt 



1 



-ir9SKK {K+K- + (-l) Is K°K°) 

JK 
Jtt 

-^(<?s w / 2 )( , V?) 2 + j^gs^d^d^T] 

J TV J 7T 

-igsKK (d„K+d»K- + {-lysd^d^K ) 

JK 

1 

7- 
1 

Ik 
1 



J 7T 



2 gs 1 KKeB fl K + <9 M K~ 



J7T 



B^B^TT+TT- 



1 

Ik 



2 gs 11 K K e 2 B ll B^K + K- 



(24) 



(interactions with rf are omitted here for brevity) . Here S 
stands for any scalar field, arj,/o or a, and P - for pseu- 
doscalar tt~ tt ,^ or K^, K°, K° and r\. We have in- 
troduced the effective couplings gs-K-n, gs-q-q, etc. listed in 
Table [2 1$ — for /o and <r and Is = 1 for ao . Cou- 
plings are expressed in terms of the model parameters Cd, 
c m and 8, see also the expression (|A2[) for the C 9jS coef- 
ficients. The Lagrangian ([24j) leads to the vertices shown 
in Fig. El 

Given this set of interaction terms, the leading con- 
tribution to the 7*7^ vertex comes from the one- loop di- 
agrams [T2J. The mechanism of the <j> meson decay via 
the kaon loop was first considered in a different formalism 
in [TU] and is consistent with the data [B]. We would like 
to stress that in the current approach the loop mechanism 
is a predicted subprocess following directly from the La- 
grangian, rather than an assumption. In particular, for the 
case of the 7r°7r°7 final state both the kaon and pion loops 
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Table 1. The vertices from Resonance Chiral Lagrangian terms (12 1 1) - (1 231) . The dashed line stands for pseudoscalar meson 
(momentum I), double solid — for vector meson, wavy line — for photon (momentum q). 



Diagramm 
Vertex function 


OA) t v \ 

=&AA> ' 

eF v [Qh'XQfi - g v ^qx) 


(M) ^ 
= « 

N 

N 

l- 

2/J LV l A V l A J 


1 + 

OA) *' O) 

N 1 
< 

l- 

Mr + z+) M + «+)x] 




p CJ cj> 


P 






P 






Up = M 
K± (fp = f K ) 
K° Up = Ik) 




2 
1 

-1 



1 
1 




-V2 
-V2 


2 
1 





1 






-y/2 





1 1 -1 

2 6 3V2 





(")< 



/+ 



2e 'g^fp 
jp 



Figure 3. The vertices corresponding to the Lagrangian (|24l) . 
The dotted line stands for a scalar meson S, the dashed one 
— for a pseudoscalar P. Couplings are shown in Tabled 

contribute. The latter are very important in the region of 
the p resonance (recall that the 7* invariant mass y/s is 
not constrained to the <f> meson mass). 

When working with the three-point vertex functions 
j*jS, we factorize the kaon-loop part in the ao case and 
separately the pion-loop and kaon-loop part for /o and <r, 
as illustrated in Fig. [U (see |Appendix C| for details). The 
7*(Q M ) 7(fc")5(p) amplitude reads 

T"" = -ie 2 (Q"k>* - gTQ ■ k)F s ^{p 2 ,Q 2 ). (25) 

The 7*(Q 2 ) — > jS(p 2 ) transition form factors (FF's) have 
the form 

Ff or ^(p\Q 2 ) = Gg. 7 (p 2 ,Q 2 ) +G^ 7 l 7 (p 2 ,Q 2 ),(26) 
iV> 2 ,Q 2 ) = G<$. y (j?,(f) + G^> 2 ,Q 2 ), (27) 



^ or > 2 ,Q 2 ) = Gif 7 l 7 b 2 ,g 2 ), 

where the terms 



(28) 



G^> 2 ,Q 
for 5 = /o, a, and 



,(Q 2 ),(29) 



«.> 2 ,Q 2 )=%^/f 44 )^;(Q 2 X30) 



27T 2 m 2 K 



K " l K 



follow from (j2"Tj) -([2"4" |) . and the pion and kaon electromag- 
netic form factors, F* m (Q 2 ) and F^ n (Q 2 ), follow from (EI]) 
and (El). The terms 



Gskk(p 2 ) = 1/flc (gsKK(m 2 K ~p 2 /2)+g SKK ) , 
Gs.Ap 2 ) = I// 2 {gs^iml - p 2 /2) + g S7T7r ) , (31) 
for S = fo,<7 and 

G ao KK(p 2 ) = (ga KK(m 2 K - p 2 /2) + g a „ifif ) , 

G ao7rr ,(p 2 ) = l// 2 (g am) (m 2 + m 2 - p 2 )/2 + <w,X 32 ) 

have the meaning of momentum-dependent SPP vertices. 
The expression for I(a, b) in (f2l?]) - (P(I|) coincides with that 
of [xTim] and for convenience is given in |Appcndix C 



The scalar meson contribution relevant to the it 
final state is 



O — _ u _u 
O, 7T 7T 



E 

S=f , <y 



D 



Ap 2 )g s ,Ap 2 ) (cg, 7 b 2 



rG^> 2 ,Q 2 







and in the tt t] case one has 

= D ao (p 2 )G ao ^(p 2 )G { a ^Jp 2 ,Q 2 ). 



(33) 



(34) 



We use the scalar meson propagator Ds{p 2 ) in the 
form E3] 

Ds\v 2 ) =P 2 - M 2 S + M s &n(f 8 , tot{M 2 ) 

+^ f Sitot (p 2 ) (35) 



with 



r to t,s{p 2 ) 



r to t,a {p 2 ) = r aa ^ v ( P 2 ) + r aa ^ KR ( P 2 ). 



(36) 



Contributions of heavy particles to the total widths, e.g., 
rf ^. vr] {p 2 ), are neglected. Modified widths r in the above 
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Table 2. Effective couplings for scalar mesons [24] (to be used 
with vertices of Fig. [3]). Model parameters are Cd and c m ; the 
scalar octet-singlet mixing angle 9 is defined in Eq. (|A4[) ; 77' 
couplings are omitted; singlet couplings Cd and c m are related 
to Cd and c m in the large- iV c approximation. Notice that the 
entries relevant to the r\ meson correct the results of Table 9 
in Ref. Q2]. 



9fvv 

9fKK 


= -2c m m 2 (2 cos 6» - \[2 sin6»)/^3, 

= -c m (2{C 2 s {2m 2 K -ml) + C 2 q ml) cosO 

+V2(C 2 (4m 2 K - 2ml) - C" 2 m 2 ) sm6)/V3, 
= -c m m 2 K (4 cos 6> + V2 sind)/V3 . 


9fvv 

9fKK 


= 2 c d {2 cos 9 - V2 sin 9)/VS, 
= c d (2(C 2 +C 2 )cos8 

— \/Z(Ly q — zo s ) smo)/ v , 
= c d (4cos6» + V2 sin 9) /V3. 


Q a 7T 7T 

QaKK 


= -2 c m ml (V2 cos 9 + 2 sin 9)/V3, 

= -c m (~V2(Cl(4m 2 K - 2ml) - C" 2 m 2 ) cos6> 

+2{C 2 s {2m 2 K - ml) + C" 2 m 2 ) sm9)/V3, 
= -c m ra|-(-^2 cos# + 4 sin6>)/^3 . 


Q a 7T 7T 
9 (J 7] 7] 

(jtrKK 


= 2c d (y2cos6» + 2 sin6»)/\/3, 
= c d (V2{C 2 -2Cl)cos9 

+2(C 2 + C 2 )sin9)/V3 , 
= c d (-V2 cos 6» + 4 sin 6»)/\/3. 


g a KK 


= — v/Scmmlf, 

= — 2^2 C q Cm ml , 


QaKK 


= V2c d , 

— 2y/2Cq Cd ■ 


9f7rr} — dftrr) — Qo-kt] — Qairrj — , 




= —igsTrir, gs-/KK = —igsKK , 

= gSnn, gSjjKK = gSKK 



expressions are defined similarly to the tree-level decay 
widths given in |Appendix B[ sec Eqs. (|B10|) . but the an- 
alytic continuation is used: 



v7¥) = e iAr3Wp))/2 V\f¥)\, 



(37) 



see Ref. [M]. 

By construction, the functions /1 in (|3"3"1) . ([M]) are of 
the chiral order 0(p 6 ): the diagrams of Fig. [4] arc 0(p A ) 
and SPP transition is 0(p 2 ). 



4 Vector contribution 

For 7* — > (• • ■ ) — > 7r°7r°7 the vector contribution mecha- 
nisms are listed in Tableland the corresponding diagrams 
are shown in Fig. [5] 

For the odd-intrinsic-parity vector- vcctor-pscudoscalar 
and vector-photon-pseudoscalar interactions we use the 
chiral Lagrangian in the vector formulation for spin-1 
fields. As shown in [57], the use of vector formulation for 



1~ fields ensures the correct behavior of Green functions 
to order 0(p 6 ), while the tensor formulation would re- 
quire additional local terms (see also discussion in the Ap- 
pendix F of [IS])- We choose Lagrangians of Ref. [2"Tll28j . 
that arc 0(p 2 ) and 0(p 3 ), for construction of the vector 
7 VP and double- vector VVP contribution to /.;. General 
Lagrangian terms are given in |Appcnclix A 



Assuming exact SU(3) case, the 7 V interaction can be 
written as 



C lV = -ef v d»B v (pl 



—ijj 
3 



^1 ^ 
— <?VJ 



(38) 



with V^y = dfjYv — <9„ V M and fy — Fy /M p is the coupling 
for the vector representation of the spin-1 fields [55] . 

The interactions of vector mesons in the odd-intrinsic- 
parity sector read 



C 



VjP 



4V2eh\ 
3/ w 



{p^ + 3^ + 3e u<j> ^)d v Tr 



+ [(3p 0fl + w")C, + 20"C»] d v r] 



(39) 



r 40-y 

*~VVP — 7 — 6 

fa 



+r][(d ft p a " 'd a p op + d^^d ^)^ C q 
V2 



,(40) 



where is the totally antisymmetric Levi-Civita ten- 

sor. As before, we omit the 77' meson. 

As it is also seen from (|3"9"]) and P0]) , the transitions 
7(/>7T°, (j)p°ir and tfiuir] are related to a small parameter 
£u4>i responsible for the uu + dd component in the phys- 
ical c/> meson. The parameter e' is responsible for the G- 
parity- violating (fiuiir vertex, caused by isospin breaking. 
The coupling constants fy, hy and 6y are model pa- 
rameters. Numerical values for all parameters are given 
in |Appcndix B| 

Due to a similar structure of the Cyp 1 and Cyyp 
interactions, the processes 7* — > VP1.2 —> P1P2I (one- 
vector-meson exchange) and 7* — > V a — > VbP\.i — > 
P1P2I (double- vector- meson exchange) can be described 
together. For this purpose it is convenient to introduce 
the form factors F 7 *yp(<5 2 ) which describe the transi- 
tions 7*(Q 2 ) — > VP including both these mechanisms. 
Of course, the vector resonance enters off-mass-shell. 
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Figure 4. Scheme for the 7*7/0 and 7*70- (top) and 7*700 (bottom) transition. Each "loop blob" corresponds to a set of 
diagrams following from the Lagrangian, as explicitly shown in [12] 



For the 7* — > Vn° transition we obtain 



For the 7r u 7r u 7 final state one obtains: 



+e UJ<j) V2<j v f v Q 2 D^Q 2 )] : (41) 
F^ n (Q 2 ) = ^-[V2h v - a v f v Q 2 D p {Q 2 ) 

JIT 

+e'^a v f v Q 2 D (f> {Q 2 )}, 
F r ^(Q 2 ) = e u ^[V2h v - avfvQ 2 D p (Q 2 )] . 



fl — "J E F 1 'Vtt{Q 2 )F 1 *Vti{^) 
V—p,uj 

x[(k-Q + l 2 )(D v (R 2 + ) + D v (R 2 _)) (44) 
+2k ■ l(D v (R 2 + ) - D V (R 2 _))], 

I2 = \ E F rVv (Q 2 )F r vA0)[Dv(Rl)+D v (R 2 _)], 



V—p,oj 

1 

4 



fs=~J E F rVn (Q 2 )F^ V7T (0)[D v (R 2 + )-D v (R 2 _)] { 



V—p,u) 



The vector meson V = p^to^cj) propagators are 

D V (Q 2 ) = [Q 2 - M 2 + iV&FtotyiQ 2 )}- 1 . (42) 
with an energy-dependent width for the p meson 



2 \ ^ V ^p 
Ftot,p(Q ) ~ 



487T/4Q2 
1 



(Q 2 -4mlf 2 6(Q 2 -4 m l 



+ -{Q 2 -Am 2 K Y' z e{Q 2 ~Am K 



(43) 



and the constant widths for the u and 4> mesons. 

In terms of these FF's we find the contribution to the 
functions /j (see Eq. 0) coming from the processes (J3J). 



where the contribution proportional to 
F 7 * ( f, 7: (Q 2 )F 1 * ff (O) oc e^. has been neglected. The 
momenta are defined as 



R% = (1/4)(Q 2 + I 2 + 2k ■ Q ± 2(k ■ I + Q ■ I)), (45) 

or equivalently R\ = (k + p\) 2 and R 2 _ = (k +P2) 2 - 
Similarly, for the 7* — > Vrj transition we obtain FF's 



F Tm {Q 2 ) = C q F ruj7T {Q 2 ) 7 (46) 

Fj^niQ 2 ) = 2 C s ^[V2h v - a v f v Q 2 D^Q 2 )} 
-s^(C v + C s )-^-a v fvQ 2 DUQ 2 )- 
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Correspondingly, the contribution to the functions /, 
for the 7r ?77 final state is 

fY = -\ E {F^vA0)F rVv (Q 2 ) 

x [{k ■ Q + l 2 )D v {R\) + 2k ■ lDv{R\)] 

-h-T-y* Vr} 

x[(k-Q + l 2 )D v {R 2 _) - 2k ■ ID V {R 2 _)] }, 



& = \ E {F~evAV)F r vn(Q 2 )Dv{Rl) 
+F 1 , v , 1 {0)F 1 , v AQ 2 )D v (R 2 _)Y 
\ E {F^vA0)F rVv (Q 2 )D v (R 2 + ) 

V— p,w,0 

-F 7 , V)) (0)F 7 .^(Q 2 ) J D y (i?2)|. (47) 



5 Numerical results 

In this section we present the numerical results obtained in 
our framework. The model-dependent ingredients, namely, 
the functions /i,2,3 are given in Sections [3] and HI 

The values of the model parameters, which we used 



in our numerical results, are listed in Appendix B The 
masses of vector and pseudoscalar mesons are taken 
from pQ . The coupling of vector mesons to a pseudoscalar 
and photon hy is estimated from the tree-level decay 
width. The scalar meson couplings and mass parameters 
were found from the fit l24l. 



5.1 Scalar mesons and <j> radiative decay 

As we discussed in this paper, in e + e~ annihilation to 
7r°7r°7 and 7r°?77 both scalar © and vector decays ([3J 
contribute to the observed events. The KLOE Collabo- 
ration has reported data on the invariant mass distri- 
butions [2|4] at y/s = M^, in which the vector meson 
contribution has been subtracted. In |24] we performed 
a combined fit of dB(tfi — > ao7 — > it°rp()/dy~p* and 
dB(4> -> (/ ,cr)7 -> it^^^/dy/p 2 to the KLOE 2002 
data 00], considering only scalar meson contributions. 
We have found the inclusion of the a meson into the frame- 
work important, and have fixed the numerical values of 
scalar meson couplings and mass parameters within the 
model, for more detail see [23]. In Fig. |H] we show our 

model results for dB(cj) S7 — >• PiP 2 -y) / dyfp 2 , eq. ([T9]) . 
at y/s = M$. In this and subsequent plots we use the no- 
tation m^o^o and m„ OT o for yfp 2 . Note that only the scalar 
meson contribution to the P1P27 final state is plotted in 
this Figure. The plot for the 7r 7r°7 final state shows a 
rather good fit g3] to the KLOE 2002 data g], where 
both /o and a are taken into account. 

In 2009 the new KLOE data [3J on the tt ?77 chan- 
nel appeared. A comparison of the model prediction for 





Figure 5. The vector, 7* — > VPi — > P1P2J, and double vector, 
7* — ?> V a — > VtPi — > P1P27, contributions 

Table 3. Mechanisms of the vector contribution. Notice that 
some of the channels, suppressed due to small parameters, can 
be enhanced in the vicinity of the corresponding resonance (e.g. 
7* — > 4> — > LJ7r°, see the text) 



Dominant 



Suppressed 



m 7 



(...) 



7r°7r°7 



1-vector 


(A U ), {"A 




2- vector 




//V u ), 
(p -><^ ) 


in 7* —►(•••) — > 7r°7;7 : 


1-vector 


( P7 r u ), (ujtv") 
(PV), M) 


(<M' J ) 

(H) 


2- vector 


(p — > ujtv"), (u) -¥ pn°) 
{p -> prj), (u -» uri) 


(p^<M u ), W>^/97r u ) 
<jyq), (<f>-¥ wq) 



4> — > ao7 — > 7r°?77 with these new data is also shown 
in Fig. [S] (bottom). We leave a refined fit of these new 
data for the future. Notice, if one adds vector contribu- 
tions to er(e + e~ — > r]TT ^) according to Table [3J then 
the shape of the invariant mass distribution, calculated 
from eq. (fl"9]l . changes: cf. Fig. [6] (bottom) and Fig. [7] 
It turns out that the 2009 KLOE data [3] are better 
described by the total contribution rather than by the 
scalar part alone. Note that in Refs. [3j[8] it was claimed 
that the <j) — > 7r°?77 decay is dominated by the <j> — > 007 
mechanism and the vector contribution is very small: 
B(e+e~ -> VP -> ?y7r°7) < 10" 6 . 



5.2 The 7* 



lutt and 7* 



ujit contribution 



For the moment, to follow KLOE analysis [5] we neglect 
the G-parity-violating vertex 4>ujtt°, i.e., we set e' = 0. For 
illustration we introduce the constant C£ w [5lll8|. This 
constant can be obtained in terms of form factors (|4"T]) 



167ra 



r(s) F r6W (0), 



(48) 



leading to 



C£ w = -167rcv 4v/ ^ y (V2h v - avfvsDp(s) 
ps (0.597- 0.542 i) GeV" 2 



(49) 



at y/s = M^. The KLOE result [5] for the same constant is 
C£ w = 0.850 GeV" 2 (y/s = M$). Thus our prediction for 
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Figure 7. Invariant mass distributions in the e + e~ annihila- 
tion to 7r°777 for yfs — M$, where the total contribution (vector 
and scalar) is taken into account (cf. Fig. [6] (bottom)). Data 
are from [3]. 
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Figure 8. Partial differential cross section of e + e~ annihila- 
tion to 7r°7r°7 for yfs = due to the 7* — > p — > ujtv mecha- 
nism compared to 7* — > (p, <j)) — > ujtt 



Rcf. 



(Eqs. (32), (33)): 



Figure 6. 

tion to 7r°7r°7 (top panel) and 7r°?77 (middle and bottom panel) 
for yfs = M$. Data are from [4] (top), [2] (middle) and [3] (bot- 
tom) 



C? v = -16™ 

J IT 



4y/2h v 
s 



(V2h x 



(50) 



i + /V 

(1.06- 0.69 i) GeV~ 2 



(D p (s) + i3 p ,D p ,(s))) 



the absolute value, \C^\ = 0.751 GeV 2 , which includes 
only the 7*(— > p) — > ujtt mechanism, is smaller than that 
of KLOE by about 15%. 



p = p(1450) 

meson which is not included in the present calculation. To 
estimate the role of the p' in the constant C£ w , we follow 



for Pp, = -0.25, M p , = 1.465 GeV, r p > (M 2 , ) = 400 MeV 
and obtain |C£J = 1.27 GcV~ 2 . 

Next we turn on the parameter e' responsible for the G- 
parity- violating (pituj vertex and check how the value 
changes. Omitting p' we have 



= -167TCV 



1 (V2h v - a v f v sD p (s) (51) 



f 2 



+ ^crvfyse'D^s)) sa (0.52 - 0.72 i) GeV" 2 
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and obtain |C£ OT | = 0.892 GcV -2 . While makingthis esti- 
mation the value s' = —0.0026 has been chosenQ. Appar- 
ently, the present model with the lowest nonet of vector 
mesons, supplemented with the G-parity-violating effect, 
allows one to obtain the value for close to the KLOE 
value 0.850 GeV -2 . Influence of the e' parameter on the 
cross section is presented in Fig. [8j 

Therefore, the difference between the C£ T value origi- 
nating from the 7* (— > p) — > ujtt mechanism, and the value 
measured by KLOE may be explained by the p' meson 
and/or G-parity-violating contribution. To clarify further 



this issue, an analysis of data at s 
sential Q. 



1 GcV^ will be cs- 
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5.3 The 7* — > <fi — > pn and 7* — > u> — > pn 
In a similar manner one can define C pn (s): 



16ira- Ry* f 



v {s) F rp7r (0) = G p7r (s) 



(52) 



where 



= -Wira 



4y/2h x 



(0.091 - 0.002 i) GeV -2 



(V2h v - avfvsDu(s)) 



(53) 



and 



C r p l s = -16™ 



4x72/11 



9/ 2 



fvs 



-0.0052. 



(54) 

The KLOE values for these constants are C r p % s ps -0.0057 

and Cp^ — 0.26 GeV -2 . However, in the experiment, 
they are entangled and one has to compare the total 
contributions. Using the values (|53p and (fM]) we have 
\C p7r (M^) \ sa 1.2, which is in a reasonable agreement with 
KLOE fit |CV(M 2 )| ps 1.3. 



5.4 Full model prediction for the cross section 

Interference of leading vector resonance contributions (pit) 
and (w7r) is presented in Fig. [9l One can see a destructive 
interference. 

The interplay of the scalar @ and vector decay (J3|) 
contributions to da/dy/jr* is shown in Fig. [TU] (for y/s = 
M^). One observes a complicated interference between 
vector and scalar contributions. We see that in the case of 
the 7r 7r°7 final state the vector contribution has the same 



cj7r deter- 



1 Of course the experimental decay width <; 
mines only the absolute value of this parameter. 

2 At s = 1 GeV 2 the G-parity-violating vertex is suppressed, 
whereas the p' mechanism survives. Therefore, any difference 
in the values of at two energies, s — 1 GeV 2 and s = Mi, 
would indicate sizeable G-parity-violating effects. 



Figure 9. Vector and double-vector decay contributions to 



da/dy 



of e 4 



7r 7r 7 at y/s = M<f, in the approximation 



£w0 = 0.058, e' — —0.0026. The (cjm) channel is negligible and 
not shown in the plot 
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Figure 10. Differential cross section dojdyff? of the e + e 
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for y/s = M4 
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annihilation to (top panel) and it t]j (bottom panel) 

for Js = 1 GeV 



size as the scalar meson one and is much smaller than the 
scalar one for the 7^777 final state. 

Notice that there exist the off-peak (*Js = 1 GeV) data 
collected by KLOE. The </> meson decays get strongly sup- 
pressed and the total cross section is determined by the 
vector contribution only. In order to support the related 
activity and provide the important model estimates, we 
include this case into our numerical calculation. The cor- 
responding results are presented in Fig. 111! 



6 Conclusions 

We presented a general framework for the model- 
independent decomposition of the differential cross sec- 
tion for the final-state radiation in the reactions e + e~ — > 
7r°7r°7 and e + e~ — > 7T 777, for which the ISR contribution 
is absent and the leading-order cross section is determined 
solely by the FSR mechanism. 

We calculated the explicit form of the functions /j, 
which carry the mo del- dependent information about the 
processes. Scalar resonance, vector and double vector me- 
son exchange contributions are considered. Notice that all 
the relative phases are fixed from the Lagrangian of Res- 
onance Chiral Theory. The only exception is the sign of 
the e' parameter, which is related to a rare <f> — > lutc decay. 



The Lagrangian is taken at the linear-in-rcsonance 
level in the even-intrinsic-parity sector and at the bilinear- 
in-resonance level in the odd-intrinsic-parity sector. For 
agreement with data, the RxT Lagrangian with the low- 
est nonet of vector and scalar mesons [5U] was extended 
by including some SU(3) symmetry breaking effects. At 
the same time, we tried to keep the number of model pa- 
rameters as small as possible, using additional constraints. 
The model parameters for the scalar sector were obtained 
from the fit [24] to the KLOE data 04]. 

As a by-product, we also obtained predictions for var- 
ious transition form factors: 7*75, SPP, 7*V~P and 
j*PP. These expressions follow directly from the La- 
grangian, and the corresponding parameters are fixed to 
a large extent. 

The numerical results for the differential cross section 
da/d\J%? are given for two cases: \fs = 1 GeV and y/s = 
Mrf, and demonstrate an interplay of the scalar and vector 
decay contributions. The influence of the scalar and vector 
contributions on the cross section is studied in detail. 

The main conclusions of the numerical studies are the 
following: 

— for the 7r°f77 final state the vector contribution is much 
smaller than the scalar one at y/s = whereas for 
the 7r°7r°7 channel the vector and scalar contributions 
are of the same size; 

— among the vector contributions to the 7r°7r°7 chan- 
nel the leading one comes from the 7*(— > (p',4>)) 
tun mechanism; comparing to the KLOE fit [3] we 
have concluded that about 85% of this contribution 
is caused by the p intermediate state, and the rest can 
be explained either by the p(1450) or by the G-parity- 
violating process: 7* — > <f> — > lot:. New experimental 
data at y/s = 1 GeV can help to clarify which of these 
two mechanisms is responsible for the rest; 

— at yfs = 1 GeV the scalar contribution is suppressed 
and the total cross section is determined only by the 
vector contribution both for the 7r 7r°7 and 7^777 chan- 
nels. 

At the end, we would like to emphasize that the de- 
veloped approach allows one to obtain the cross section 
and branching fraction close to the experimental results. 
The main advantage of this approach is a small number 
of model parameters. 

The proposed framework can be implemented in a 
Monte Carlo generator, for the inspection of the com- 
pletely differential characteristics of the reaction, and thus 
is useful for a data analysis and a detailed comparison of 
various models. 
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Appendix A: Pseudoscalar mesons, scalar 
multiplet and the RxT Lagrangian 

In chiral theory, the pseudoscalar mesons n, K, 77 can 
be treated as pseudo-Nambu-Goldstone bosons of spon- 
taneous G = SU(3) L x SU(3) R to H = SU(3) V broken 
symmetry. The physical states 77, 77' can be introduced 
using the scheme with two mixing angles (9q, $s); for a 
review see [30]. The adopted scheme is consistent with 
chiral theory and takes into account the effects of U(l) 
axial anomaly and SU(3) flavor breaking (m s 3> m u ,d)- 
In our notation [23] the pseudoscalar nonet reads 



Interaction terms for the pseudoscalar and vector 
mesons read 



(Al) 



/ ■K°+c qV +c" q n' 



exp 



where 



V2U 



7T+ \ 

V2 _ /if 

'-K- -C.r, + CWJ 



c n = 



c s 
c. 



y/3 cos (9$ — 


0o) 


u 




y3 COS ($8 — 


0o) 


u 




v / 3cos(0s — 


Oo) 


u 




v / 3cos(0g — 


&o) 



■ cos 



h 



sm ( 



(A2) 



— V2cos# 8 + — sin# 

JS JO 

i V2 cos O + 7- sin 8 

JO J8 

i cos 6» 8 - — -v^sin^o 
h Jo 



The vielbein field which represents the pseudoscalar 
: iu + D p u + and \+ = u + x u+ + 



mesons is u,, 



u\u is the explicit symmetry-breaking term, x ~ 
diag(m^, m^., 2m|- — rri^) in the isospin symmetry limit. 



f 

T v 



2^2 



(WO 



(V^'V), (A3) 



here (• • • ) stands for the trace in flavor space. 

For scalar mesons we assume the nonet symmetry of 
the interaction terms and multiplet decomposition 



' ao =Sz, 
fo =5o cos 9 — Ss sin#, 
(7 =Sq sin 9 + Sg cos 9, 



(A4) 



where S3 is the neutral isospin-one, Sg is the isospin-zero 
member of the flavor octet. The angle 9 is the octet-singlet 
mixing parameter, and a = ./o(600). The interaction La- 
grangian for scalars takes the form 



source, F, 



1"' 



d v B p is the electromagnetic field 
tensor. It appears in the chiral covariant derivative, which 
in our case is reduced to D^U = d p U + ieB p [U, Q] and in 
the tensor = eF^ v (uQu + + u + Qu), where the quark 
charge matrix Q = diag(2/3, —1/3, —1/3). 

For calculations in the even-intrinsic-parity sector, 
the leading-order RxT Lagrangian for pseudoscalar, 
scalar, vector mesons and photons was derived by 
Ecker et al. [2D]. The spin-1 mesons are described by an- 
tisymmetric matrix tensor fields V 1 ^ and this Lagrangian 
is equivalent to the ChPT Lagrangian at order 0{p A ) 
(see [20,25 for details). In our application we have some- 
what released the rigor of RxT and use different masses 
of resonances (M p ^ M u ^ M<p and M a ^ M ao ^ M /o ) 
without specifying a pattern of flavor symmetry breaking 
(cf. Ref. EI]). 



^scalar = C d (Su^) + C m (Sx+) 



(A5) 



There are known problems with a rigorous inclusion of a 
and /o(980) into any RxT multiplet [3T]. However, there 
is also a number of successful applications [321I33] of a 
scheme similar to (|A4[) . In studies of <j> radiative decays 
this scheme was also applied in [2~4"ll26j . 

Due to nonet symmetry, the relation for scalar sin- 
glet So and octet S oct coupling constants holds, c mt dS = 
c m , d (S oct + So/V3) ■ In nomenclature of Ref. [2D] this 
relation implies c m d = c m .d/V3- 

In the odd-intrinsic-parity sector the flavor SU(3) sym- 
metric Lagrangian [271128] reads 

C odd = h v e p „ afs (v»{u v ff + ffu v )) (A6) 
Appendix B: Model parameters 



The electromagnetic field B^ is included as an external Masses 



The following values for the meson masses arc used 
in our numerical calculations [T]: M p = 775.49 MeV, 
M u = 782.65 MeV, M4 = 1019.456 MeV, m n = m w ± = 
139.57 MeV, m,o = 134.98 MeV, m K = 493.68 MeV, 
m v = 547.75 MeV. 



Mixing parameters 

The values of the 77 mixing angles 9q = —9.2° ± 1.7° and 
6» 8 = -21.2° ±1.6° are used [35], thus / 8 = (1.26±0.04)/ 7r 
and f = (1.17 ± 0.03)/ w , where f„ ~ 92.4 MeV. Thus, 
one obtains C q ps 0.738 and C s ps 0.535. 

The uj(j) mixing is given by one parameter e W = 
0.058 [3D] . The states of "ideal mixing" = (uu+dd) / y/2 
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and <pid = ss arc expressed in terms of the physical ones 
(mass eigenstates) as 



Uid — OJ + £u<f><P, 
4>id = 4> - £lj<I>U. 



(Bl) 



In order to include a G-parity-violating ifiujir vertex we 
determine the parameter e' from the <f> — > urn decay width: 



ri 



1 9<t>- 



' TT I P-jz 



12tt 



(B2) 



where P7I = y/(M$ + m\ - M*) 3 /(4 Mj) - m*, the effec- 
tive coupling in our formalism is = Aaye' / f v . Using 
the experimental value for the <f> — > ujtt decay branching 
ratio B = (4.4±0.6) x 10" 5 [Mj and a v = 0.34 one obtains 
le'l = 0.0026. 



Couplings in the even-intrinsic-parity sector 

The condition Fy = 2 Gy for the model couplings is 
used in our calculation to make the one-loop amplitude 
finite jT2] without use of counter-terms. This relation has 
been addressed in |25| in a different context, namely it 
has been shown that the constraints imposed by the high- 
energy behavior of the vector and axial-vector FF's lead 
to it, in addition to the relation FyGy = f 2 . Note that 
Fy = 2 Gy also appears in alternative models, e.g., Hid- 
den Local Gauge Symmetry Model and massive Yang- 
Mills theory for vector mesons, see a discussion in [251 ■ F° r 
numerical calculations we use Gy = f,r/V2 = 65.34 MeV, 
F v = 2 G v = 130.68 MeV. 

Alternatively, respecting phenomenology, one may fix 
Fy and Gy by means of fitting the measured partial decay 
widths of the vector mesons (see, e.g., [20]) at tree level. 
In particular, for p — > e + e~ one has 



leads to h 



V 



0.041 ±0.003. 



One can use a special short-distance constraint of 
R%T in order to relate ay to fy and hy. Namely, one 



can require the form factors 



to vanish at Q 



—00 

as expected from QCD. In this connection we refer to [37l 
138] . where in the framework of RxT high-energy behav- 
ior of three-point Green functions VVP, VAP, AAP has 
been studied. 

At Q 2 -> -00 the propagators (g2J) D V (Q 2 ) -> 1/Q 2 
and we obtain the following relation (neglecting mixing) 



\2hy — cryfy = 0. 



(B6) 



This constraint reduces the number of independent pa- 
rameters in the model, in particular, expresses the poorly 
known parameter ay via hy and fy, which can be fixed 
from data. Thus we obtain ay ~ 0.34. 

Notice, an additional constraint on the parameters ay, 
hy and fy follows from the short-distance behavior of the 
7*7*7T° form factor (see a discussion in Ref. [55]): 



N c 1- 

^ + mV2hyfy 



8&vfv 



(B7) 



It allows to further reduce the number of independent pa- 
rameters. For example, one can leave fy to be the only 
independent parameter and deduce from (|B6[) and (|B7p 



Nr. 



ay 



32 TT 2 f 2 



32y/2 TT 2 / V 



(B8) 



r. 



V 



p—ye-i 



Yl-KMn 



(B3) 



and for the p —> itit the tree level width is given by 



r, 



r 2 

1 2 



487T/4 



Ami 



(B4) 



The experimental data are the following [T]: r(p 
TT+7T-) = 146.2±0.7MeVandr, 



p— 



7.04±0.06keV. 

Values obtained in this way are Gy = 65.14 ± 0.16 MeV 
and Fy = 156.41 ± 0.67 MeV. The estimated values sup- 
port the Fy w 2 Gy conjecture. 



which results in the numerical values ay = 0.329 and 
hy = 0.0395 — fairly close to those obtained with the use 
of Eq. (TBBl . 

In favor of broken flavor SU(3) symmetry, one may 
introduce separate couplings for each vector meson, i.e. 
replace fy by f p , f 
h 



and further hy by h lpTr , h 1U}Vl 



(7 VP transition) , and also ay by cr wpT , a ppv , 



Parameters for scalar mesons 



Couplings in the odd-intrinsic-parity sector 

The coupling constant fy is given by fy — Fy / M p sa . 1 7 . 



The parameter hy can be fixed from the V 
width, in particular, the p — > ttj width 



r(p -> 7T7) 



AaM 3 p h v 
27/2 



Mj 



Pj decay 
(B5) 



The widths for ao — > 77 and fo — > 77 decays are expressed 
in terms of ([2^]) - ([3H]) . for example: 



4 4 

e p 
QAtt^Jp 2 



\G a K u P 2 M 



(B9) 
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Table 4. Scalar meson parameters [24]. Couplings and mass 
parameters are given in MeV. 



the former with the contact 7* — > tt + tt coupling and the 
latter with an intermediate vector resonance. They read 



Cd Cm 


M ao M fo M a 


6> 


93ti l 46l y 2 


1150^ 986.1^;* 5041^ 


36° ± 2° 



The strong decay widths of the scalar mesons in the lowest 
order (tree level) are 

r / 2\ _ I G aQ7r ri (p 2 ) 1 2 



^^rTtloop 



7 — yV— >7rloop 



4e2 * 1 z\ 2 t 

n G fomr (p )—I 



(4tt) 2 1 " J0 ""^ ' m a 

xD fo (p 2 )G fo7ITr (p 2 ) 

~4e 2 z ^ 2 2 

(47r) 2 r i Wo-lPJ TO 2 

xD fo (p 2 )G fo7ITT (p 2 ) 

x ^F V G V Q 2 D P (Q 2 ). 



f P 2 



(C2) 



? 2 p 2 



(C3) 



j(p 2 +m 2 — m 2 ) 2 
4p2 



x^p 2 /4-m 2 , 



ZL! m 2 ; (BIO) The o meson propagator D r (Q 2 ) is given in (|42|) . The form 

factor G/ o7r7r (p 2 ) is given by (|3T|). The loop integral I(a, b) 
can be found, e.g., in [TT] and [H], and reads 



87rp 2 



where p 2 is the invariant mass squared of the scalar meson; 
see also definition (|31[) , ((32)) . For discussion of momentum- 
dependent couplings Gspp{p 2 ) and constant SPP cou- 
plings of other models (e.g., [TU]) see Ref. [2"5] . 

The finite-width effects for scalar resonances are very 
important and expressions (|B9|) . (|B10|) do not have phys- 
ical meaning of decay width, when evaluated at the res- 
onance peak value of p 2 . Nevertheless, in several papers, 
e -g-> [T21l2"Ull3"5] . the tree-level expressions of a similar form 
were used to find the model parameters (cd, c m and 9) 
from measured widths. It was observed [2~4"ll26[[4"0] that 
the coupling constants could be better determined from 
fitting the 7T7t and Ttr/ invariant mass distributions in 
e + e~ — > <j) — > 77T7T and e + e _ — > 4> — > 77r?7 reactions. The 
fit results [53] are shown in Tableland these values are 
used in our numerical calculations. Notice that for this fit 
we used data from [3] (7r°7r°7) and [5] (tt°t]j). Recently, 
a new KLOE result for the latter appeared [3], and we 
find reasonable agreement with it without refitting, see a 
discussion in Section [5] 



I{a,b) 



2(a-b) (a -by 



(a -by 



9[ b 



9 ~ 

a 



(C4) 



with 



arcsm 



(2^) 



x > 



log 



i+(x) 



in 

2 



, < x < |, 
x < 0, 



g(x) = { \^T=Ax~ 



\/4x — 1 arcsin^7^= 

1 „ "+(*) . 



X > 



.(*) 



, < s < \, 
x < 0, 



n±(aO = (l ± VT 3 ^) 

ZX 



(C5) 



For a reference, we remind the alternative notation of |12) : 
>f(m 2 ,p 2 ,Q 2 ) = (a-b)I(a,b), 



l/(Q-k) = 2/(Q 2 -p 2 ), 

with a = Q 2 /to 2 and 6 = p 2 /m 2 . 

Combining (|C2|) and (|C3[) . one finds 



(C6) 



Appendix C : Example of the factorization of 

the 7* — y 7/0 — > 77r°7r° transition amplitude 

In this Appendix we sketch the general structure of the 
scalar meson contribution ff giving emphasis on the ap- 
pearance of the electromagnetic form factors of the pseu- 
doscalars in the formulae ([55]) and (|34[) . 

Consider the part of the M^ v amplitude of 7* — > 
7/0 - > 77r°7r° with a pion loop transition, . Fig- 

ure S] is of help and one observes two terms 



TT loop 7— IK loop 



•7T loop' 



(CI) 



A/7 

TZLOOp 



-4e 2 i -o. 2 /Q 2 p 2 

r l ^/oTTTrl-P Jo I ™2 ' ™2 



/ ^ \9. 1 JOTTTi V-T / 9 1 9 5 9 

(47r)^ m4 \ to w m w 

x^ /o ( P 2 )G /oT7r (p 2 )F; m (g 2 ) 

EEE- le 2 rr^/ b 2 )G /07r7r b 2 )GW, 7 (p 2 ,g 2 ), 

where the two-photon form factor of a scalar meson 
G^» 7 (p 2 , Q 2 ) is given in ([2"5]) . The pion electromagnetic 
form factor F^ m {Q 2 ) in RxT is given by 



F? m (Q 2 ) = l-^Q 2 D p (Q 2 ) 

J 7T 



(C7) 
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■:'+S 



V 



Figure 12. The 0(p ) electromagnetic vertex of a (off-mass- 
shell) pseudoscalar meson in RxT . All possible intermediate 
vector resonances V = p°, ui, </>, ... in general contribute. For 
real photons only the first term on the r.h.s. is non-zero. 



Factorization in the part of the M^ v amplitude {7} 
with a kaon loop transition, M^ lo , is analogous. The 
kaon form factor in RxT is 



1 - (\D P m 

6 3 



(C8) 



The vector meson V = p,u,(f> propagators are given 
by (gll). The form factors in form (fC7]l and (|C8]) include 
contributions from the photon-vector transition (vector 
meson dominance, VMD) and the direct jPP interac- 
tion, sec Fig. [T21 The detailed discussion of two versions 
of VMD (VMD1 and VMD2) is given in the review gj. It 
turns out that the RxT corresponds to the VMD1 version. 

For discussion of the one-loop modification of the elec- 
tromagnetic vertex and RxT -motivated calculation of the 
kaon form factor see [12] . 
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